Legolas 2.0: Improvements and extensions to an MHD spectroscopic framework

We report on recent extensions and improvements to the Legolas code, which is an open-source, finite element-based numerical framework solve linearised (magneto)hydrodynamic equations for a three-dimensional force- thermally balanced state with nontrivial one-dimensional variation. The standard Fourier modes imposed give rise complex, generalised non-Hermitian eigenvalue problem solved quantify all linear wave of given system in either Cartesian or cylindrical geometries. now supports subsystems eight MHD equations, allowing pure hydrodynamic setups, only density/temperature/velocity variations, option treat specific closure relations. discuss optimisations internal datastructure solvers, showing considerable performance increase both execution time memory usage. Additionally code has capability fully visualise eigenfunctions associated multiple dimensions, we apply Kelvin-Helmholtz Rayleigh-Taylor instabilities hydrodynamics, thereby providing convincing links between stability analysis onset non-linear phenomena.